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Abstract 

The present theories of galaxy formation are reviewed. The relation between peculiar 
velocities, temperature fluctuations of the microwave background and the correlation 
function of galaxies point to the possibility that galaxies do not form uniformly 
everywhere. The velocity data provide strong constraints on the theories even in the 
case when light does not follow mass of the universe. 

1. Initial Conditions 

The universe contains a wide dynamic range of objects : from stars (l M@) 
all the way to superclusters (10 16 Mq). A major question that we are unable to 
answer yet is whether the formation of structure has started with smaller masses 
clustering on ever larger scales *), or whether extremely large structures formed first, 
then subsequently fragmented into smaller ones 2 ). If we knew the precise initial 
conditions then the present structure of the universe could be derived by applying 
the laws of physics. Let us summarize, what has to be known about the initial 
conditions for this ambitious project. 

The fluctuations are likely to be adiabatic, since the specific entropy of the 
universe, ub / ti 1 is tied to microscopic parameters of particle physics. Entropy fluc- 
tuations, once popular, can be generated by huge amounts of shear, e.g. In the infla- 
tionary theories quantum fluctuations arise in a natural way. However, the necessary 
amplitude seems to require rather special prescriptions for the effective potential 3 ). 
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The initial perturbations are expected to be scale free, therefore their Fourier 
amplitude depending on the wavenumber k can be well described by a power law, 
|6jc| 2 a k n . If the spectral index is n = 1, the amplitude of the different pertur- 
bations is the same when their wavelength equals the horizon size. This ’double 
scale-invariant’ is called the Zeldovich spectrum, and is known to arise in inflation- 
ary scenarios 4 ^. 

There are severe constraints on the fluctuation amplitudes. If the fluctua- 
tions were adiabatic, the perturbations of the metrics generate fluctuations in the 
temperature of the microwave background. On small angular scales (4.5 arc mins) 
these limits are extremely small 5 ):AT/r < 2.9 x 10~ 5 . The H-He plasma becomes 
gravitationally unstable only after recombination, at Z ~ 1000. At this point the 
density and temperature fluctuations are similar, 3 A T/T = A p/ p. Since the stan- 
dard growth of fluctuations in a flat universe is (1 + Z) -1 , this does not leave enough 
margin for fluctuation growth, the fluctuations cannot reach the nonlinear stage our 
universe seems to be in today. Present calculations confirm 6 ) that if the universe is 
baryon dominated, only prohibitively high initial fluctuation amplitudes can result 
in the formation of galaxies. If the universe is dominated by some form of collision- 
less dark matter, the dark matter fluctuations are unaffected by pressure, therefore 
grow even before recombination. After recombination these curvature perturbations 
caused by the dark matter will accelerate fluctuation growth in the baryons, so the 
AT /T constraints are less stringent. 

Though the initial spectrum is a power law, by the time it becomes nonlinear 
it will be considerably modified. When the universe is radiation-dominated, fluctu- 
ations within the horizon have a minimal increase 7 ), whereas the ones outside the 
horizon grow. This effect will bend the slope of the spectrum from n to n — 4 for 
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wavenumbers higher than fc eg , corresponding to the size of the horizon when the mat- 
ter and radiation energy densities were equal. The presence of the collisionless dark 
matter results in distortions of a different kind: the free motion of particles erases 
structures smaller than the free streaming scale 8 ’ 9 ’ 10 ’ 11 ). The mass scale of this col- 
lisionless damping process can be expressed in terms of the mass and entropy of the 
particles the dark matter consists of. M x « 2.2 m 8 ra“ 2 . In the case of neutrinos this 
mass takes the value of M um = 3.2 x 10 15 m 3 0 2 M©, corresponding to the comoving 
length scale A„ m = 41 m^ 1 Mpc. Depending on what the ’temperature’ of the dark 
matter is, this damping scale can change from the above 41 Mpc to extremely small 
values. The neutrinos are hot particles, since their average momentum is close to 
that of the background radiation photons. Most other candidates for the dark mat- 
ter like axions and photinos - yet undiscovered - would have decoupled much before 
the neutrinos, having a lower entropy or temperature, so they are called cold . They 
hardly move at all, their damping scale is negligible. Intermediate candidates, like a 
gravitino of 1 keV mass would be warm. 

A major underlying assumption in calculating most consequences of a given 
fluctuation spectrum is that the phases of the individual Fourier components are 
random, ie. the perturbations are a random Gaussian process. One can envisage sce- 
narios, where this will not be the case, like perturbations originating from strings 12 ). 
For a given spectrum combined with the assumption of random phases one can cal- 
culate the distribution of mass fluctuations, density of local peaks, density profiles 
around local peaks, the distribution of peaks of a given size, etc. 

The expansion of the universe is characterized by three quantities: H = pf Pcriu 
the density parameter, Ho, the Hubble constant, Ao, the cosmological constant. If 
A 0 = 0 and H = 1 the universe is flat, which appears to be necessary for inflation. 
Aq is generally assumed to be negligable. Calculations of the primordial 4 He and 
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D+ 3 He abundance indicate 13 ), that the baryon density of the universe at the time of 
primordial nucleosynthesis lies in the range of 0.01 < fl B < 0.1. This suggests that 
if baryons dominate the mass density then the universe is open by a large margin. 

Fluctuation growth also depends on the density of the universe. If fi < 1, 
the growth of perturbations effectively stops at the redshift Z = ft -1 . The detailed 
predictions of AT/T are just below the current limits if the dark matter consists 
of neutrinos with about 30 eV mass, and restrict 0 if the cold particles dominate 
the universe 6 ): 11 > 0.2 x h" 4 / 3 where H 0 = lOOh km /s Mpc. In deriving this 
limit it was assumed that galaxies follow the mass distribution: the amplitude of the 
fluctuations today was normalized to J3, the integral of the galaxy-galaxy correlation 
function £ g (r). 

2. Nonlinear structure 

Here we would like discuss the expected structure of the universe if the dark 
matter is either hot, warm or cold. Once the first mass scale in a spectrum with a 
large damping cutoff (hot) reaches nonlinearity, particle trajectories cease expanding 
away from each other and converge, resulting in the temporary formation of caustics. 
The density becomes very high and a flat ’pancake’ is formed 2 ). At first they arise at 
isolated spots where the initial velocity perturbations had the largest gradient. Soon 
these regions grow, turning into huge surfaces which intersect, forming the walls of 
a cell-structure which is itself gravitationally unstable. The methods of catastrophe 
theory were applied 14 ) to analyze structure that develops in such potential motion. It 
was found that the two dimensional pancakes are only the lowest order singularities; 
other singular topological structures should also appear. String-like features are one 
example, and they can be seen in the N-body simulations. 



When the intersection of trajectories takes place, gas pressure builds up, the 
velocity of the collapsing gas exceeds the sound speed and a shock wave is formed 2 ). 
The gas is shock-heated up to keV temperatures and cools by emitting radiation over 
a broad spectrum. The UV and soft X-ray emission can photoionize the intergalactic 
medium, making galaxy formation in regions that have not yet formed pancakes more 
difficult, which would accentuate the contrast in galaxy density between the strings 
and pancakes vs. voids, even though the density contrast may be only 3-10. 

If the dark matter is cold, then the mass autocorrelations are logarithmically 
divergent towards the smallest scales. These objects will collapse first, the scales 
determined by the baryon Jeans mass at around recombination. Collapse of larger 
scale systems follows subsequently. It is believed that the statistical properties on a 
given mass scale can be reasonably well understood by studying the Gaussian ran- 
dom fluctuations obtained by filtering out all the smaller scale contributions from the 
power spectrum. Recently as major effort has been undertaken 15 ), where galaxies 
were associated with peaks of a given height of the random fluctuation field and var- 
ious properties like correlation functions, mean shapes, densities etc were calculated 
in a manner similar to previous work on pancakes 16 ). If the dark matter is warm, 
it will still form pancakes, though of galactic size. There the cooling is much more 
efficient 17 ), those timescales will determine the fate of each object. 

In either of the above scenarios it seems to be very hard to avoid strong initial 
explosions and rapidly cooling shocks, which compress the gas and provide seeds 
for the next generation of explosions, as suggested by Ostriker and Cowie 18 ). The 
complicated nature of such calculations has yet prevented a very detailed discussion, 
but the importance of these processes is unquestionable. 
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3. Peculiar velocities 

If we knew all the parameters listed above, it would be relatively easy to 
follow the evolution of the universe- Only gravitational forces act on collisonless dark 
matter so one can numerically solve the transport equations, even in the nonlinear 
regime. This has indeed been done, as we discuss here. Given the initial conditions, 
these numerical experiments can tell us the mass distribution in the universe. One 
can hope, that the structure obtained this way will resemble the real universe, ie. 
galaxies trace the mass distribution. 

Starting from the above mentioned initial conditions extensive N-body simu- 
lations 19 » 2 °) were made. The free parameters of the calculations are H, H 0 and the 
initial amplitude of the fluctuations. For a given H one can use conservative limits 
for the age of the universe to obtain a value of H 0 . If fi = 1 , then t 0 > 12 Gy 
requires Ho < 54 km/s Mpc. The initial amplitude can be defined in various ways. 
For simulations with hot dark matter the epoch of galaxy formation Zqf was the 
redshift when 1 percent of all particles have gone through a ’caustic’. For cold dark 
matter, due to the growth of nonlinearity, £(r) is rapidly increasing both in slope and 
amplitude, just like for hot dark matter. One can define today when the correlation 
function of the particles most resembles that of the galaxies, ie. a power law with a 
slope -1.8. 

£(r) = (r/r 0 ) -1 ’ 8 

The simulations have encountered a major difficulty : the random velocity dispersion 
of galaxies is well known 21 ): 

< v\ 2 > 1/2 ?s 300 - 400km/s. 
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Both in the neutrino and cold dark matter simulations, when the density correlations 
are just about right, velocity dispersions are in the 1200 km/s range, clearly too high. 

< v 2 12 > 1 / 2 « (1200 km/s) n°' e £o 

where £o is the value of the mass autocorrelations at 5 Mpc radius. Comparing this to 
the data, this suggest that ft < 1, forbidden by the A T/T constraints. Since a low ft 
model is ruled out, the only remaining possibility is to have £o = |<5p/p| 2 fairly small. 
Then we are in a sharp contradiction with the observed galaxy autocorrelation. 

Here one should note, though, that all calculations so far have assumed, that 
the distribution of galaxies follows the mass distribution, ie. £ g (r) = £ m (r). It is 
£ m (r), which determines both the A T/T fluctuations and the peculiar velocities, and 
it is £ g (r) that we can observe. Since £ g (r) seems to be too large to be in agreement 
with either AT /T or (vJ 2 )j and changing ft does not resolve the problem, the next 
possible solution may be that the mass fluctuations are relatively small, whereas 
£ ff (r) ^ £m(0* This means, that galaxies do not form with uniform probability 
everywhere, the formation rate is ’biased’ towards some regions. 

This can be quite natural, though, since galaxies consist mostly of baryonic gas 
capable of emitting and absorbing radiation. These dissipative processes, strongly 
density and temperature dependent, occur at a different rate at different places 17 ). 
All these effects, combined with possible shock waves due to the finite pressure in 
the H-He gas, may have an important role in determining where galaxies form. As a 
result, the galaxies may not follow the light at all, so the mass autocorrelation should 
not be compared to the galaxy autocorrelation. Galaxy formation, as long as it is a 
random process, initiated by gravitational infall will be likely to start at the regions 
of highest densities. One can therefore associate the particles in these regions with 
galaxies. This ’biasing’ of galaxy formation towards these high densities is a heuristic 
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procedure, but probably a fair approximation to what really happens. The physical 
explanation of what the threshold of the selection should be is much less clear, it can, 
only be adjusted to the observed number density of galaxies. This ’biasing’ process 
enhances the correlations, without invoking large peculiar velocities. 


If we consider the large scale velocity fields, they provide strong upper limits 
to the ’biasing’ factor 22 ). The dispersion of the center-of-mass velocity of a sphere 
with radius R is given by 

{V 2 ) = (Ho/) 2 [°° dk\6 k \ 2 W 2 (kR) 

Jo 

where / — fi 0 * 6 and W(kR) is the window function, the Fourier transform of the 
spherical distribution 23 K The window function effectively eliminates contributions 
from scales smaller than R , so { V 2 ) is a genuine measure of the large scale fluctua- 
tions, which are believed to be still close to linear. During the last few years there 
were several attempts to determine the peculiar velocities of spheres of galaxies cen- 
tered around us, although the errors are considerable 24 ). The results are not yet 
conclusive, but potentially they are an important test of the fluctuation amplitude. 

Another measure of the large scale structure of the universe is the cluster - 
cluster correlation function £ cc (r). It has the same functional form as the galaxy 
autocorrelation, but the amplitude is considerably larger 2S ); 


£cc(r) = ( 


36Mpc 


- 1.8 


Furthermore, the amplitude is dependent upon the richness class. It has been shown 
recently, that this richness dependence can be nicely explained, if we assume that the 
universe has a scale invariant property over the volume of the Abell catalogue. For 
each cluster sample one can derive the mean distance between clusters (D = n -1 / 3 ), 



which would uniquely characterize the richness. If we measure the length in these 
units, the richness dependence disappears 26 ) : 

£ cc (r) = 0.35(r/I>)- 18 

Originally Mandelbrot 27 1 has suggested such a ’fractal’ structure for the universe. 
The physical meaning of this scale invariance is not absolutely clear. It is unlikely that 
it could be generated via nonlinear gravitational dynamics, since the corresponding 
velocities on 40 Mpc scales would be enormous. There are suggestions, that cosmic 
strings may have such an effect 28 ), but other explanations attribute the difference in 
the clustering amplitude to the fraction of galaxies associated with clusters 29 ). 

Recent calculations indicate, that for certain kinds of fluctuation spectra the 
correlation function of ’biased’ regions may be a power law over a wide dynamic 
range, and the slope of the power law would depend on the threshold set for galaxy 
formation 30 ), contrary to previous work, claiming that £(r) would be amplified by 
a constant factor 31 ). 

5. Conclusion 

All the present theories of galaxy formation fail to explain the observed uni- 
verse in its full complication. The recent observations of the microwave background 
fluctuations provide the strongest constraints on present theories. The details of 
galaxy correlation properties are a new challenge, indicating that galaxies are un- 
likely to be tracers of the mass distribution. The peculiar velocity field of galaxies 
and clusters may provide a way to probe the fluctuations even in this case. 

I would like acknowledge useful discussions with Simon White, Dick Bond, 
Jim Bardeen and Lars Jensen. 
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